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ABSTRACT 

We calculate the evolution of 6 Li generated from cosmic rays produced by an early population of massive stars. 
The computation is performed in the framework of hierarchical structure formation and is based on cosmic star 
formation histories constrained to reproduce the observed star formation rate at redshift z < 6, the observed 
chemical abundances in damped Lyman alpha absorbers and in the intergalactic medium, and to allow for an 
early reionization of the Universe at z ~ 1 1 by Pop III stars as indicated by the third year results released by 
WMAP We show that the pregalactic production of the 6 Li isotope in the IGM via these Pop III stars can 
account for the 6 Li plateau observed in metal poor halo stars without additional over-production of 7 Li. Our 
results depend on the efficiency of cosmic rays to propagate out of minihalos and the fraction of supernovae 
energy deposited in cosmic rays. We also compute the cosmic ray heating of the IGM gas. In general, we find 
somewhat high temperatures (of order 10 5 K) implying that the cosmic rays production of 6 Li may be required 
to be confined to the so-called warm-hot IGM. 

Subject headings: Cosmology - Cosmic rays - Big Bang Nucleosynthesis - Stars: abundances 



1. INTRODUCTION 

The Lithium observed in low metallicity environments such 
as the atmospheres of halo stars offers a unique probe into 
two very distinct mechanisms of nucleosynthesis: the big 
bang and cosmic rays. Big bang nucleosynthesis (BBN) 
produces predominantly 7 Li while cosmic ray nucleosynthe- 
sis (CRN) produces roughly equal numbers of the 6 Li and 
'Li isotopes. The standard lore tells us that the bulk of 
the Pop II 7 Li abundance is produced by BBN with addi- 
tional contributions at the 10% level being supplied by galac- 
tic CRN (or GCRN). BBN lays down t he primordial abun- 
dance which dominates the Spite plateau ( Spite & Spite 1982) 
while GCRN supplies a metallicity dependent supplement 
which would result in a small slope in Li vs. Fe regressions 
jRvan etal. l l2000). Li on the other hand, produced only in 
GCRN would be expected to show a strong (log-linear) re- 
lation with Fe, with an abundance which is a fraction of the 
7 Li ab undance. Indeed the first observat i ons of 6 Li at r Fe/H] 
2 dSmithj Lambert. & Nissenlll993t _ | Hobbs & Thorburn 
" " Smith Lambert. & Nissen 1 119981 ICavrel etaL 
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1999; Nissen et al. 2000) wer e at t he expected level of 
6 Li/ 7 Li ~ 0.05 (ISteigman et al. lll993t iFields & Qliv"e1ll999T: 
Vangioni-Fla m et al. 11999ft . 

This simple picture has been shaken by rece nt observations 
leading to tw o distinct Li problems. WMAP (Spergel et al. 
120031 12006ft has determined the baryon density of the Uni- 
verse to high accuracy, Og/i 2 = 0.0224 ± 0.0009, corre- 
sponding to a baryon-to-photon ratio, 77 = (6. 14 ±0.25) x 
10~ 10 . At this value of 77, th e BBN predicted value of 7 Li 
" x 10" 10 (iCvburt. Fields. & Olive ll200!l 
, 7 Li/H = 4.9t\i x 10 " 10 llCuoco et al. I 



is 7 Li/H = 4.27+1;^ 
20031 ICvburt l f2004) 
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2004), or 7 Li/H = 4.15^ x 10" 10 (ICoc et al. I 12001 . 
These values are all significantly larger than most deter- 
minations of the 7 Li abundance in Po p II stars which are 
in the range 1 - 2 x 10" 10 (see e.g. ISpite & Spitel 119821 
iRvan et a l. 112000) ) and even larger t han a recent determina- 
tion by Melendez & Ramirez (2004) based on a higher tem- 
perature scale. The second Li problem originates from re- 
cent observations of 6 Li ( Ascl und et all 120061) which indi- 
cate a value of [ 6 Li] = log 6 Li/H+ 12 = 0.8 which is inde- 
pendent of metallicity in sharp contrast to what is expected 
from GCRN models. In effect, we are faced with explain- 
ing a 6 Li plateau at a level of about 1000 times that ex- 
pected from BBN (iThomas. Schramm. Olive. & Fie ldsll 19931 
Vangioni-Flam et al. I1999I) . Here, we will concentrate on the 
latter of the two Li problems. 

Different scenarios have been discussed to explain 
the abundance of 6 Li in metal-poor halo stars (MPHS). 
ISuzuki & Inoue I J2002) discussed the possibility of cos- 
mic rays produced in shocks during the formation of the 
Galaxy ; which wa s consistent with 6 Li data available at that 
time. iJedamzikl ( ^000| l2004albl). Kawasaki et alJ J2005I) 
Jedamzik et al. (2005), Kusakabe. Kaiino & Mathews 
(2006), and Pospelov (2006) consider the decay of relic 
particles, during the epoch of the big bang nucleosynthe- 
sis, that can yield to a larg e primordial abundance of 6 Li. 
IFields &"Pr odanovic ( 2005 ) have studied in detail the lithium 
production in connection to gamma rays, using a formalism 
similar to ours. 

Previously, we computed the evolution of the 6 Li abun- 
dance produced by an initial burst of cosm ological cosmic 
rays (CCRs) (Rollinde. Van gioni. & Olive 1 120051) rhereafter 
RVOI]. We found that the pregalactic production of the 6 Li 
isotope can account for the 6 Li plateau observed in metal poor 
halo stars without additional over-production of 7 Li. The de- 
rived relation between the amplitude of the CCR energy spec- 
tra and the redshift of the initial CCR production put con- 
straints on the physics and history of the objects, such as Pop 
III stars, responsible for these early cosmic rays. Due to the 



2 



subsequent evolution of 6 Li in the Galaxy through GCRN, we 
argued that halo stars with metallicities between [Fe/H] = -2 
and -1, must be somewhat depleted in 6 Li. 

Here, we will employ detailed models of cosmic chemical 
evolution to derive the total CCR energy as a function of red- 
shift. In parti c ular, w e will make use of models discussed in 
iDaiene et al. I (|2006) which were constructed to account for 
the observed star formation rate (up to z ~ 6), the observed 
SN rates (up to z ~ 1 .5), early reionization of the intergalactic 
medium (IGM) (at z ~ 11), and reasonable chemical abun- 
dances in the interstellar medium (ISM) and IGM of proto- 
galaxies. 

These models allow us to determine the total energy in- 
jected in cosmic rays as a function of redshift. Once a spectral 
shape for the source function is given, we must account for the 
propagation of CCRs. Furthermore, we must distinguish be- 
tween diffusion in the ISM and propagation in an expanding 
IGM. We first consider the resulting nucleosynthesis if most 
of the high energy particles are ejected from the minihalos 
where star formation begins and escape into the IGM. This 
will lead to a physical picture similar to the assumptions made 
in RVOI. In this context, we will see that cosmic -ray heating 
of the IGM may place important constraints on the scenario. 
Depending on the fraction of SN energy deposited in CRs and 
the efficiency of shock acceleration of CRs leading to their 
escape from structures, it may be necessary to consider in ad- 
dition to the IGM production of 6 Li, the in situ production 
when CRs are confined to structures. 

Our paper is organized as follows: In the next section, we 
wi ll briefly describe the cosm ic chemical evolution models 
of IDaiene et al. I f2004l 120 06). In section 3, we describe the 
general features of cosmological cosmic rays related to the 
production of Lithium in the IGM. There we also present our 
detailed calculations of the 6 Li abundance based on the Pop III 
production of CCRs in the IGM, at high redshift. In section 
4, we comment on the in situ production of 6 Li in Pop III 
minihalos. We summarize the status of the CCR origin of 6 Li 
in section 5. 

2. SN HISTORY AND COSMOLOGICAL COSMIC RAYS 

The cosmic star formation histories considered are based 
on the detailed mo dels of chemical evolution derived in 
Daign eet al. I J2p06). The models are described by a bimodal 
birthrate function of the form 

B(m,t,Z) = M*")Mt)+M'")M z ) (!) 
where 4>i(2) is the initial mass function (IMF) of the normal 
(massive) component of star formation and i/>U2) is the respec- 
tive star formation rate (SFR). Z is the metallicity. The normal 
mode contains stars with mass between 0.1 M Q and 100 M 
and has a SFR which peaks at z ~ 3. The massive component 
dominates at high redshift. The IMF of both modes is taken 
to be a power law with a near Salpeter slope so that, 

(j)i(m) oc (2) 

with x = 1.3. Each IMF is normalized independently by 

r a' sup 

/ dm m(pi(m) = 1 , (3) 

Jm m[ 

differing only in the specific mass range of each model. Both 
the normal and massive components contribute to the chem- 
ical enrichment of galaxy forming structures and the IGM, 
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FIG. 1. — The star formation rate for Model 1, and the rapid burst Model 
le. The dotted curve shows the SFR of the massive mode of Model 1. Data 
come from Hopkins i2004). 



though the normal mode is not sufficie nt for acco unting for 
the early reionization of the IGM (Daigne et al. 2004). 

Here, we r estrict our att e ntion to the best fit hierarchi- 
cal model in Daigne et al. (2006) in which the minimum 
mass for star formation is 10 7 M . The distribution of 
struct ure masses is based on the Press- Schechter formalism 
JPress & Schechterl[l974llJenkins et al. 11200 11) . The normal 
mode SFR is given by 

tpi (f ) = viMstrm exp {-t/n ) , (4) 

where M struct is the mass of the structure (which includes dark 
matter), t\ =2.8 Gyr is a characteristic timescale and v\ = 0.2 
Gyr -1 governs the efficiency of the star formation. In contrast, 
the massive mode SFR is defined by 

4>i{t) = v z Mi m exp (-Z IGM / Zen, ) , (5 ) 

where Mism is the mass of the baryonic gas in the structure, 
with h>2 maximized to achieve early reionization without the 
overproduction of metals or the over-consumption of gas. We 
adopt Z cr it/Z Q = 10~ 4 . 

We also restrict our attention to Model 1 of Daigne et al. 
(2006) to describe the massive mode. In Model 1, the IMF 
is defined for stars with masses, 40 M Q < m < 100 M R . 
All of these stars die in core collapse supernovae leaving a 
black hole remnant. The coefficient of star formation, V2 is 80 
Gyr -1 . Star formation begins at very high redshift (z ~ 30) but 
peaks at a redshift z — 12. Note that the absolute value of the 
SFR depe nds not only on vj, b ut also on the efficiency of out- 
flow. See Daigne et al. (2006) for details. We also consider 
an example of a model (with an IMF as in Model 1) in which 
the massive mode occurs as a rapid burst at z — 16 designated 
as Model le. Thus, most Pop III SN occur within a very short 
period of time at this redshift. This model is therefore sim- 
ilar to the SFR assumed in RVOI. In Figure ^ we show the 
SFR, for Models 1 and le (including the normal mode). Also 
shown by the dotted curve is the SFR for the massive mode 
alone in Model 1 . 

We do not he re consider mod els with very massive stars in 
Pop III. In D aigne et al. l f2006). two very massive mode mod- 
els were considered in addition to Model 1 described above. 
In Model 2a, the massive mode consists of stars with masses 
in the range 140-260 M Q , and in Model 2b, the range used was 
270-500 Mr. In Model 2b, there is significant metal enrich- 
ment llHeger & Wooslev 12 002) which results in a diminished 
SFR and hence fewer cosmic rays available for 6 Li produc- 
tion. In Model 2b, one expects total collapse and little or no 
production of cosmic rays. 
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The rate of core collapse supernovae (SNR) can be calcu- 
lated directly in terms of the IMF and SFR 



SNR = 



x(SM Q , 



dm<f)(m)^p{t — t(ot) ) , 



(6) 



ininit) 



where m mjn ^ is the minimum mass of a star with lifetime, r, 
less than t. When a star undergoes core collapse, the mass 
of the remnant is determined by the mass of the progenitor. 
We assume that all stars of mass m > 8 M e , will die as su- 
pernovae. For stars of mass 8 M© < m < 30 M Q , the remnant 
after core collapse will be a neutron star ofmwl .5 M . Stars 
with 30 M < m < 100 M become black holes with mass 
approximately that of the star's helium core before collapse 
(Hege r et al. 1120031) . We take the mass of the Helium core to 
be 

M He = ^-(m-2OM ) (7) 



for a star with main sequence mass m (iHeeer & Wooslevl 
2002). The supernova rate ultimately determines the metal 
enrichment of the ISM and when coupled with the model of 
outflow also determines the metal enrichment of the IGM. 

The energy emitted in each core collapse, E cc corresponds 
to the change in gravitational energy, 99% of which is emitted 
as neutrinos. In the cases where collapse results in a neu- 
tron star, E cc = 3 x 10 53 ergs. For stars that collapse to black 
holes, E cc is proportional to the mass of the black hole. For 
masses less than 100 M , we take E cc = 0.3Mn e ■ We will 
parametrize the energy injected in cosmic rays per supernova 
as 

eE cc (m) 



£? C r('«) 



100 



(8) 



where e is the fraction of energy in the remaining 1% (i.e. en- 
ergy not in neutrinos) deposited into cosmic rays. Given the 
IMF described above, the massive mode is dominated by 40 
M stars for which the total energy per SN in CRs is §m = 
10 52 ' 8 e/// ergs. In contrast, the normal mode associated with 
Pop II, is dominated by lower mass stars for which the energy 
per SN in CRs is Sn = 10 515 e// ergs. While it is quite plausi- 
ble that e// and e/// differ (indeed we would expect e// < eui), 
we will for simplicity assume en = £/// = e = 0.01 -0.30 as a 
broad and conservative range. 

The SNR derived from Eq. [6] is shown for both Models 1 
and le in Fig. |2] (lower panel). In the upper panel of Fig. [2] 
we show the energy density in CRs injected per year. The CR 
production rate in Model le (shown by the dot-dashed curves) 
is similar to that Model 1 (shown by the solid curves) below 
a redshift of about 10, as would be expected from the SFRs 
shown in Fig.[Q Dashed lines show the rates of CRs generated 
by massive Pop III SNe alone, while dotted lines correspond 
to Pop II SNe ejection. The energy density in cosmic rays is 
dominated at large redshift by Pop III SNe due to the corre- 
sponding IMFs and mass range associated with the two modes 
and the dependence of E cc on the progenitor mass. 

The metallicity evolutio n in both th e IGM and the ISM has 
been derived bv iDaigne et al. I J2006I) and is shown in Fig. [3] 
For Model le, the metallicity in the ISM rises very quickly 
to [Fe/H] ~ -2.5, whereas in Model 1, the initial enrichment 
occurs rapidly only to [Fe/H] ~ -4. The IGM abundance are 
several thousand times smaller. Both models have the same 
metallicity as a function of redshift below z ~ 15 in the ISM. 

3. PRODUCTION OF LITHIUM IN THE IGM 
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FIG. 2.— CR history predicted bv lDaigne etaTI |2"0(%1) . The SNR (lower 
panel) and energy density in cosmic rays (upper panel) are shown in the case 
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In RVOI, we considered a single burst of CRs whose total 
energy was fixed in an ad hoc way so as to reproduce the ob- 
served 6 Li abundance. We now relate the production of CRs 
to the detailed model of cosmic chemical evolution in which 
the SN history is completely determined by the SFR and IMF 
of the model and are constrained to reionize the IGM at a 
redshift z — 11, match the observed SFR at z < 6, as well as 
chemical abundances at z < 3. As a consequence, the energy 
density in cosmic rays is determined by the model and we 
can derive the abundance of Lithium produced in the IGM as 
described below. 

3.1. Formalism of Cosmological Cosmic Rays 
3.1.1. Ejection ofCCR into the IGM 
The total kinetic energy given initially to CRs by SN is 

^sn(z) = (1+z) 3 / dmcj>{m)i>(t-T{m))£'cR(m). 



max(SMQ ,rn min ^ 

(9) 

where it is understood that the appropriate IMF is used for 
computing CR energy density due to Pop II or Pop III SN. 
For example, using Eq. [9] we can estimate from Fig. ^ the 
Pop III contribution to the CR energy density. In Model 1, 
if we take ip — 0.02 M Q yr" 1 Mpc"^ from z = 10-30, and 
approximate <j> by a delta function at m = 40 M , we obtain 
<?sn(10) ^ 7.8 x 10" 13 e ergs/cm 3 . Note that the result of the 
full calculation of Eq. [9] is a factor of about 5 larger than 
this, because the contributions of more massive stars (more 
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massive than 40 Mq) can not be neglected. This corresponds 
to <?sn — 3.6 x 10 e ergs/cm 3 . For Model le, we see that the 
burst is very intense, with a SFR which reaches 20 M yr -1 
Mpc~ 3 , though the duration is only about 3 x 10 6 yr. In this 
case, we approximate that <£sn(16) — 1.6 x 10 _11 e ergs/cm 3 . 
The full calculation yields #sn — 2.5 x 10 _11 e ergs/cm 3 . 



3.1.2. Source Spectrum 

We next briefly describe our treatment of cosmological cos- 
mic rays, their transport and the mechanism for the produc- 
tion of 6 Li, in the case of one burst of CRs at a given red- 
shift z n - Our form alism is directly derived from the work of 
iMontmerlel d!977l) . hereafter Mil. We begin by defining the 
source function for our CR distribution. The source function 
Q(E,z) is defined as a power law in momentum, 

Q{E,z) = C(z ) - ^-p f^tGeV-WV 1 ], 
(£(£ + 2£ )) (7 V 

(10) 

where C(z) is fixed by our normalization of the source func- 
tion to the total energy from both Pop II and Pop III SNe using 



^SN(z)=^SN,n(z) + <*SN,m(z)= / EQ(E,z)dE (11) 

10 6 GeV 3 . We set 



min -= 0.01 MeV and E max 



We take £, 
7 = 3 in this section. 

The efficiency of the SN to eject CRs outside the structure 
and into the IGM depends on many different physical param- 
eters. Roughly, low energy particles will lose all their energy 
inside the structure and in this section, we take a simplified 
approach where the CRs spectrum is simply cut at a given 
energy, E cut and decreased by a constant factor e s hock- 



Qvm{E) = 



if£<£ C ut 

^shock _ ^CUt 



(12) 



In lDaigne et al. I J2006I) . the efficiency of the baryon outflow 
rate coming from the structures is dependent on the redshift. 
It accounts for the increasing escape velocity o f the structure 
as the galaxy assembly is in progress. In fact, in Daigne et al. 
(2004) two sources of outflow were considered : a global 
outflow powered by stellar explosions (galactic winds) and 
an outflow corresponding to stellar supernova ejecta that are 
pushed directly out of the structures as chimneys. However, 
velocities in ISM gas are of the order of a few 100 km/s while 
CRs are mostly relativistic. Thus, the CR ejection processes 
considered here should be independent of the overall outflow 
of gas and heavy elements. As we will see, the production 
of Lithium is proportional to the energy available, and thus 
to e s hock- Most of the CR production will occur at large red- 
shift, when Pop III are dominant (especially in the case of 
the Model le). For simplicity, we will assume that all CRs 
are ejected from structures and adopt a constant value for 
^shock = 1- We comment on the possibility of e s hock < 1 in 
§4. 

3.1.3. Production of Lithium in the IGM 



3 Note that in RVOI we used £ m ; n = 10 MeV. When the integral diverges, 
i.e. when the power law 7 is larger than 2, this can affect the result by a factor 
as large as 2. The cut-off of 0.01 MeV is related to a minimum kinetic energy 
for CRs to escape the star itself. 



If Ni(E,z) is the comoving number density per (GeV/n) of 
a given species at a given time or redshift, and energy, we 
define Nj.n(E,z) = N(E,z)/nn(z), the abundance by num- 
ber with respect to the ambient gaseous hydrogen (in units 
of (Gev/n) -1 ). The evolution of Ni.n is defined through the 
transport function 



dN 'V 9 (UA7 N N W n 

-^r + dE {bNi '« )+ i»= Qi >«- 



(13) 



where 7b is the lifetime against destruction and b describes 
the energy losses due to expansion or ionization processes 
((Gev/n) s" 1 ). The energy and time dependencies can be sep- 
arated as b(E,z) = -B(E)f(z). We can distinguish two cases 
depending on whether losses are dominated by expansion or 
by ionization. The general form for the redshift dependence, 
whe n expansion dominates is fp.(z ) = (1 +z)~ 1 |dz/d?|7/Q 1 
(e.g. lWick. Dermer & Atovanll2004l) . Other contributions to 
B or /, do not depend on the assumed cosmology and are 
given explicitly in M77. 

In M77, two important quantities, z*(E,E',z) and E^(E,z) 
are used in this formalism. Given a particle (a or lithium) 
with an energy £ at a redshift z, z*(E,E',z) corresponds to 
the redshift at which this particle had an energy E'. E s (E,z) 
is the initial energy required if this particle was produced at 
the redshift of the burst, z s . In particular, z*(E,E s ,z) = z s - 
The equation that defines z* (Eq. A5, M77) is 6z*/dE = 

-[B(E)f(z)\dz/dt\]- 1 (dz*/dz). 

The evolution of the CCR a energy spectrum is derived, 
using Eq. A8 of M77 



4> a {E) 13 4> a {E s ) 
n° 0' <l> a {E) 



exp (-0 



1 



\b(E, Zs )\ \dz*/dE'\ E , 
(14) 

where & a p(E,z) = & a {E ,z)/nn(z) is the flux of a's per co- 
moving volume 



<S> a p(E,z)=PN a p(E,z) 



(15) 



and (3 (/?') is the velocity corresponding to energy E (E s ); £ 
accounts for the destruction term (Eq. A9, Mil). The CR 
injection spectrum, 4> a is proportional to the source spectrum, 
Q- 

The evolution of CR flux and the production of Lithium can 
be computed step by step directly with the transport function 
(Ea.ll3>.The abundance by number of lithium (/ = 6 Li or 7 Li) 
of energy E, produced at a given redshift z, is computed from 



dN LH (E,z) 
dt 



a aa ^i (E,E')n He (z)$ a fl(E',z)dE' 



= a,(E)K ap ^ a (4E,z) [(Gev/n)" 1 s" 1 ] , (16) 

where a aa ^i(E,E') = ai(E)S(E-E'/4) and K ap = 0.08 is 
the abundance by number of 4 He/H. We use cross sections 
based on recent measurements related to the a + a reaction 
and provide a new fit for the production of 6 Li and 7 Li 
jMercer et al. "l200l . 

If we now consider the contribution of each individual burst 
at each redshift : 

(Li/H)( Z ) = Li/HBBN+J^dz'jf dz" 

x S E a,(E)K ap ® a (4E,z")dE \dt/dz"\ \dt/dz'\ 

(17) 
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3.2. Results 

Although the exact evolu tion of CR confinement is difficult 
to estimate, Ensslin (2003) discusses th e relation of the diffu- 
sion coefficient with a magnetic field. Jubelgas et al. (2006) 
propose a simple model where this coefficient varies as the 
inverse of the square root of the density. The efficiency of the 
diffusion will then decrease with the density (as could be in- 
tuitively inferred). The n, at large r edshift, the structures are 
smaller, less dense (e.g. Zhao et al. 2003) and the primordial 
magnetic field could be expected to confine less than it does 
today, which corresponds to a large value for e s hock- As noted 
above, we will assume that c^q^ — 1, bearing in mind that 
this approximation should not be valid at low redshifts (z < 3). 
In effect, our results for the production of Li will depend on 
the product of the CR escape efficiency, e s h ck, and the effi- 
ciency for converting SN energy into CRs, e. We will return 
to discuss the value of e s h oc k further below. 

In this work, we have also introduced the energy cut-off pa- 
rameter E cut . Its influence is actually straightforward. If E cut 
is below 10 MeV, it does not modify our results for the Li 
abundance. Given the shape of the Li production cross sec- 
tion, only a particles around 10 MeV will produce Lithium. 
If particles have a lower energy initially, they will never be 
available to create Lithium. While if they are at higher energy, 
they will loose their energy during their transport through the 
medium, and at some time reach the optimum energy. Thus, 
if the cut-off is larger than 10 MeV, there will be a delay until 
some particles loose enough energy to produce Lithium effi- 
ciently. As shown in RVOI, the production of Li decreases 
very rapidly with time due to the global expansion, as a result 
this delay is unfavorable to Lithium production. 

Pop III production of CRs dominate over Pop II at high red- 
shift and most of the production of 6 Li is due to the Pop III 
SN. Model le is closer to our RVOI assumption of a single 
burst. In Fig.0] we show the evolution of the Li abundances 
for e s hock = 1 in Model 1 and Model le as labeled. The CR 
efficieny, e has been fixed at e = 0.04 so that in Model le, 
assuming that E cut < 10 MeV, the total amount of Lithium 
produced at z = 3 is [ 6 Li] =-1 1.2. This value is perfectly con- 
sistent with the observations of 6 Li in halo stars. In Model 
le, the sum of all energy within the burst of CRs at z ~ 16 
corresponds to 10~ 12 ergs cm -3 (for e = 0.04), and according 
to RVOI (6.3 x 10" 13 at z = 10) would result in producing 6 Li 
at a value similar to that in the plateau. 

In Model 1 , the total energy, integrated over the full star for- 
mation history, is only about a factor of 3 less than that found 
in Model le. Consequently, the production of Lithium is re- 
duced to e( 6 Li) = -1 1 .7 at z = 3, for e = 0.04 and thus is con- 
sistent with observations for an increased value of e = 0.15. 
These values are consistent with expectations that roughly 
10% of non-neu trino SN energy is converted to CR accelera- 
tion JPrurv et al.lll989l) . For e = 0.15, the energy density in 
cosmic rays in Model 1 is 5.4 x 10~ 13 ergs cm -3 which is still 
slightly lower than the estimate in RVOI for a burst at z = 30 
and very similar to the one needed at z = 10 to produce the 
6 Li plateau. Note that, as claimed in RVOI, 7 Li is not over- 
produced by this process. 

In Fig.|4] we also show the Pop III contribution alone to the 
Li production (dashed curves) for both Models 1 and le. As 
one can see, Li production is largely dominated by the Pop III 
contribution, though at lower redshifts, the production from 
the normal mode is non-negligible. If our calculation is ex- 



tended to z = 0, we obtain a small enhancement in the 6 Li 
abundance as shown by the thin curves in Fig. 0] However, 
as noted earlier, below z — 3, we expect that as the structures 
are larger and contain more baryons, CR escape will be lim- 
ited resulting in a smaller value for e s h oc k- In this case, we 
expect the in situ production of Li in the ISM to dominate as 
discussed below in §4. 

For E cut = 100 MeV, the production would be decreased by 
one order of magnitude. Had we chosen a spectral index 7 = 
2, our results for both Models 1 and le would be diminished 
by a factor of about 40. 
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FIG. 4. — The production of Lithium in the IGM by CCRs assuming E cllt < 
10 MeV and 

e shock — 1 as a function of redshift for both isotopes in both 
Models 1 and le. Here, e is fixed at 0.04 so that Model le produces log 6 Li/H 
= -11.2 at z = 3. Since the Li abundance scales with e, adopting e = 0.15, 
would increase the Model 1 abundance to -11.2 as well. The contribution 
from Pop III stars alone to the Li abundance is shown by the dashed curves. 

The injection of CRs in the IGM will heat the medium. 
In fact, CCRs were predicted to heat the IGM and 
thus avoid the problem of overcooling i n the IGM gas 
JBlanchard. Valls-Gabaud. & Mamonlll992l) Following the 
analysis of Samui, Subramanian, & Srianand 1 (120051) . we find 
that the temperature reaches as high as 10 7 K in Model 1 at 
z ~ 5 when E cut < 0.1 MeV and 

e shock — 1 -0. There is a rela- 
tively strong correlation between the induced temperature and 
the CR energy cut-off as seen in Fig. [5] where the temperature 
due to CR heating is shown as a function of redshift for three 
choices of E cut . For E cut = 10 MeV, we see that the temper- 
ature is held to the range 10 4,5 - 10 5 K at z = 0. Because 
we have fixed e = 0.04 in Model le (as opposed to 0.15 in 
Model 1), we find somewhat lower temperatures for Model 
le as seen in Fig. [(5] In addition, the lower temperature is par- 
tially due to the difference in the SN history: in Model le, 
there is a sudden heating of the IGM that can cool for a longer 
time than in Model 1 for which heating is more progressive. 
Observations of absorption lines in quasar spectra at z ~ 4 set 
a conservativ e upper limit on the temperature of the IGM 
of 10 s K fe.g.lSchave et al.l2000tlRollinde. Petitiean. Pichonl 
l200UlTheuns et al. 120021) . Thus both Models 1 and le, with 
£cut = 10 MeV and e sn ock = 1-0 will not overheat the IGM. 
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FIG. 5. — The induced temperature by CCR heating in the IGM for Model 
1 for three choices of the CR cut-off energy, E cut as indicated. 
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FIG. 6. — As in Fig.|5]for Model le. 



be confined to the warm-hot IGM (WHIM). Correspondingly, 
the production of 6 Li in the IGM may also occur in the WHIM 
embedding the structures. In fact, the temperatures shown in 
Fig.Ela re quite r epresentative of the temperatures found in the 
WHIM JCen & Ostrikerll999tlsImcoe et alJ2002i) . The heat- 
in g of cluster gas by CRs used to produce 6 Li was considered 



al. (2C 



3.3. Discussion 



A key question pertains to the propagation of CRs into 
the IGM and the degree to which the Lithium produced is 
accreted onto the Galaxy. We assumed that all CRs with 
energies above E cut escape the structures and proto-clusters. 
Hence we have a homogeneous flux of cosmic rays in the 
IGM. More likely, in an inhomogeneous model, our flux 
would be contained in the WHIM. This scenario allows the 
Lithium produced to be accreted later, during the formation of 
the Galaxy. It then provides a prompt initial enrichment (PIE) 
at z ~ 3 required to explain the observed abundances. We can 
add this PIE to the standard GCR production of Lithium as 
in RVOI. In Fig. we show the resulting evolution of 6 Li 
as a function of [Fe/H]. The upper curve (solid blue) corre- 
sponds to Model 1 with e = 0.15 and does an excellent job 
of fitting the observed 6 Li abundances at low metallicity. The 
corresponding curve for Model le is nearly identical when e 
is taken to be 0.04. A standard GCR model of 6 Li production 
without a Population III enrichment is shown in Fig. [7] for 
comparison. 
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The IGM temperatures shown in Fig.[3]and Fig.|6]should be 
viewed as an upper limit to the temperature in the IGM pro- 
duced by CR heating. Firstly, to produce Fig. [5] and Fig. [6j 
we assumed e s h oc k = 1.0 at all times, including z < 3. At low 
redshift, the Press-Schechter formalism breaks down as struc- 
tures are no longer representative objects. As larger galax- 
ies form, the ability to eject CRs diminishes and we expect 
^shock to decrease. For e s h oc k = below z = 3, the tempera- 
tures would be lower by a factor of about 2 at z = 0. 

Secondly, our computation assumes that CRs propagate 
into the IGM in a homogeneous way. This is certainly not 
the case, as we expect the cosmic ray density and heating to 



FIG. 7. — The evolution of Li/H vs. [Fe/H] for standard galactic cosmic 
rays alone (shown by the solid green curve) and the addition of a PIE by the 
production of Lithium in the IGM as in Model 1 (shown by the solid blue 
curve) with e = 0.15. Model le would yield a nearly identical result with 
e = 0.04. 6 Li data come from AsDlund et al. 1 2006 ) with the exception of 
the lowest metallicity star which comes from Inoue et afl 120051) . 

As in RVOI, we require some depletion of 6 Li at higher 
metallicities. This should not be particularly surprising since 
all models GCRN predict a linear growth of [ 6 Li] versus 
[Fe/H]. In standard GCRN, the 6 Li abundance matched the 
observations only for [Fe/H] ~ -2. This fact was often asso- 
ciated with an energetics problem concerning the production 
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of 6 Li (Ramatv et al Jl 1 9971 l200(I iFields et aTll200l . Pop III 
production of 6 Li can successfully account for the 6 Li plateau 
at low metallicity and together with GCR production accounts 
for the present-day meteoritic value. As a result, we have 
overproduced 6 Li at [Fe/H] ~ — 1. This could be accounted 
for by depletion which is expected to be important at higher 
metallicities. It is also plausible that the degree of depletion 
increases with increasing metallicity (iPiau 12005b . 

The large temperatures produced in the IGM, prohibit the 
homogenous production of 6 Li (if£ cut < lOMeV). If one con- 
siders that the 6 Li production processes occur in the WHIM, 
the constraint on the temperature are very much relaxed. In- 
deed, this medium is denser, located around the galaxies or in 
filame nts of clusters of gala x ies. Hydrodynami c simulations 
(e.g. ICen & Ostriker 1 119991 iDave et al. Il200"!t iKan g et al. I 
2005) have long since shown that the WHIM, heated by 
shock waves during structure formation, displays two differ- 
ent phases: a cold phase, at T< 10 s K and a hot phase, at 
T~ 10 5 " 7 K. The cold one is observed through the O VI ab- 
sorption lines whose width distribution constrains the temper- 
ature (e.g. Ber geron et al. ll2002l) . The hot one may be re- 
sponsibl e for the soft X-ray background, observed by Chan- 
dra (e.g. lNicastro et al. 120031 iMc Kernan et al. 120 03) and its 
temperature can be constrai ned by the observat ion of broad 
H I lines (e.g. with FUSE. iRichter et al. 1120041) or of O VII 
lines with Chandra (Fang e^ al. 120021) and, in the future, with 
XEUS (e.g. IViel et alTI 12003). Consequently, the tempera- 
tures derived in our scenario, when taken in the context of the 
WHIM, are acceptable. In addition, in a dense region, other 
processes, such as H2 line coolin g, could low er the temper- 
ature even more. Note that in ISoltan I J2006I) hydrodynami- 
cal simulations indicate that a substantial fraction of baryons 
in the universe remains in a diffuse component WHIM. This 
component is predicted to be at T^ 10 5_7 K, as noted previ- 
ously. 

Finally, one should bare in mind that we set e s hock equal to 
a constant value (1.0). The exact value of e s hock depends on 
the strength of the magnetic field at each epoch, on the den- 
sity of the ISM etc... It is difficult to accurately estimate the 
escape fraction, which is certainly dependent o n the energy of 
the pa rticles (roughly modeled here by £ C ut)- Ensslin et al. 
(2006) have begun an investigation into the propagation of 
CRs in the ISM taking into account all physical components. 
To date, they are only able to obtain the mean values of the 
physical parameters, but this is certainly a path to follow. 
Note that if the efficiency is much less than one, the induced 
IGM temperature is reduced, as is the initial enrichment of 
6 Li. This will be partially offset because the density in the 
WHIM is larger than the mean density of the universe (used 
in the above formalism), and it is quite plausible that CRs are 
trapped there and increase the initial enrichment of Lithium. 

4. PRODUCTION OF LITHIUM IN THE ISM 

The model discussed in the previous section was a homoge- 
neous model in which the IGM flux of CCRs was controlled 
by the SN history and in particular, our choices of e//, e/// and 
as discussed above, e s hock- If CR propagation into the WHIM 
is efficient, i.e. e s hock = L their interaction with the medium 
will produce the required amount of Lithium as a PIE, as dis- 
cussed in the previous section. 

In fact, we expect CR propagation to be limited by diffu- 
sion, particularly at lower redshift so that they are concen- 
trated in or near the structures leading to not only enhanced 



6 Li production but also heating of the eventual intra-cluster 
medium as well as the warm-hot IGM. In this section, we will 
outline the computation of 6 Li production in the structures of 
the hierarchical formation scenario that end with our Galaxy. 
We assume that the CR energy output of SN is constant 
over a sufficiently long time so that we can work in the 
context of a leaky box model. We can then adopt the physical 
formalism develope d for LiBeB production in our Galaxy 
(Meneguzzi. Audouze & Reeves 1971 ), hereafter MAR, 
( Vangioni-Flam et al. I 120001 IFields. Olive. & Schramm! 
119941) . 

In GCRN, it is common to normalize the flux of CRs by re- 
producing the observed abundance of Beryllium at present. In 
contrast, here, we rely on the same cosmological models de- 
scribed above when we consider the cosmological evolution 
of CR and the 6 Li abundance. In addition, we must take into 
account the evolution with redshift of several physical param- 
eters defined below. 

Finally, the production of Lithium in the IGM as described 
in the previous section acted as an effective prompt initial en- 
richment of 6 Li in the IGM which by subsequent accretion 
and growth of structure led to the 6 Li plateau observed in our 
Galaxy. In the case of ISM production, each structure along 
the hierarchical tree inherits the medium as modified by SN 
that explode in the past. This explains the relation between 
metallicity and redshift seen in Fig. [3] Thus, when we ob- 
serve Lithium in a star at a given metallicity, we must use the 
Lithium abundance present when this same metallicity was 
reached. 

4.1. Formalism 

As before we compute the total rate of energy density put 
into CRs from SN, (?sn, using Eq. |9] We use the same source 
term, Q(E,z), defined in Eq. |K)| In the framework of a diffu- 
sion model in a medium of density p, it is useful to consider 
the scaled source q = Q/p. The normalization of Q is related 
to the energy density inside the structures. However, (?sn cor- 
responds to the energy density provided by SN ;/ it was uni- 
formly distributed within the universe, which was true when 
considering the diffusion into the IGM previously. Since this 
energy is now confined within the structure, the density there 
must be larger by a factor pj fbPb, where pt, is the average 
baryon energy density in the Universe and ft, is the fraction 
of baryons found in structures. Therefore source function is 
normalized using 



Eq(E,z)dE 



^sn(z) 

fb(z)p c 



(18) 



Note that all parameters that vary with redshi ft are give n with 
the hierarchical model provided by Daigne et al. ( 2006). 

By interaction with the particles present inside the medium, 
CR particles lose their energy. The rate of energy lost 
is noted b. We upda t e the relations in MAR using 
Mannheim & Schickeiser ( 1994). Note that this rate depends 
on the ionization fraction of the medium, x e . We use here 
x e = 0.01 but we checked that one could go up to x e = 0.1 
without modifying our results. It is also convenient to define 
the quantity w = b/pv where v is the velocity of the CR parti- 
cle. Because of the physical properties of the medium and the 
presence of a magnetic field, CRs will be confined to some 
characteristic escape length, A which in GCRN may range 
from 10 - 1000 gem -2 . Under these condition, the solution of 
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R(E')-R(E) 

A(z) 



(19) 



the diffusion equation (MAR) for the flux <J> a is 
* a (E,z) - ±fdE'q a (E,z)exp(- 

The term q a = K Qp g (evaluated appropriately in the ISM) is 

the source of a particles; and R(E) = / dE/w is the ioniza- 
tion range which characterizes the average amount of material 
that an a particle with energy E can travel before ionization 
losses will stop it. g is a function that corresponds to the dif- 
fusion solution. 
The rate of Lithium production of Lithium is then 



(20) 



Thus, if the physical properties of the structures (A,x e ) as well 
as those of the SN ejection processes (7) do not evolve with 
time, 

w)w=w(w)r 7 ^ r ?«b (2D 

H/ J z f b {z)p c Sl b dz 

The result of the integral is shown in Fig. [8] for both models 
considered here. 
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FIG. 8. — The cumulative ratio of the energy deposit ed in CRs to the mass 
of the structure as determined by the integral in Eq. 12 1 1 (curves are as in 
Fig.HJ. 



4.2. Discussion 

The main parameters of this scenario are A and 7. A is 
likely to be typically 10 gem -2 in our Galaxy. At larger red- 
shift, as discussed in 33.21 A is expected to be smaller due to 
the evolution of both the density of the structure and of the 
amplitude of the magnetic field. A lower value of A implies 
that CRs escape faster out of the structure, and thus produce 
less Lithium. This can be seen directly in Eq. \l9\ Indeed, 
one expects a strong correlation between A and the parame- 
ter e s hock discussed earlier with very small values of A corre- 
sponding to values of e s hock ~ 1 ■ 

The evolution of the abundance of 6 Li can computed di- 
rectly as a function of metallicity (see above). As discussed 
earlier, this scenario, based on the Press-Schechter formalism 
will likely break down at low redshift. If one follows the pro- 
duction of 6 Li up to that point z ~ 3, one can in fact place a 
limit to A to avoid the over-production of 6 Li in structures. In 
Model 1, we find that A < 10~ 5 g cm -2 . The upper limit in 
Model le is a factor of 4 times larger. These values are so low 
that our choice of e s hock = 1 appears to have been well justi- 
fied. Subsequently, we expect A to increase as the structure 
evolves into a galaxy where standard GCR becomes impor- 
tant. 



5. SUMMARY 

The observation of a 6 Li plateau in halo stars at low metal- 
licity poses a challenge to standard models of 6 Li production 
via Galactic cosmic ray nucleosynthesis. The level of the 
plateau is about 1000 times larger than the standard big bang 
nucleosynthesis value and about a factor of 10 larger than the 
GCRN value at [Fe/H] = -3. 

In RVOI, we showed that an early burst of cosmic rays in- 
jected into the IGM would produce a prompt initial enrich- 
ment of 6 Li. There, we used the observed plateau value to 
normalize the energy density of CRs. Here we applied this 
mechanism to a detailed model of cosmic chemical evolu- 
tion. The model was designed to reproduce the observed star 
formation rate at redshift z < 6, the observed chemical abun- 
dances in damped Lyman alpha absorbers and in the inter- 
galactic medium, and to allow for an early reionization of the 
Universe at z ~ 1 1 as indicated by the thir d year results re - 
leased by WMAR As a consequence, Dai gne et al. I (120061) . 
was able to compute the supernova rate as a function of red- 
shift. This SNR was employed here to compute the resulting 
energy density in cosmic rays. 

Our results depend on the efficiency to which cosmic rays 
are accelerated out of the first star forming structures. We 
found that for effici encies, eshock = L the Models 1 and le 
discussed in Daigne et al. (2006) can produce a plateau in 6 Li 
at the right abundance level at low redshift when the SN en- 
ergy deposited in CRs is about 4-15% of the available kinetic 
energy. This con forms well with models of CR acceleration 
( lDrurvetal.1119891) . As in standard GCRN, 6 Li depletion at 
higher metallicities is necessary if the plateau is found to per- 
sist. By computing the CR heating of the IGM, we conclude 
that most of the CR propagation and hence 6 Li production 
should be confined to the warm-hot IGM. 

We also compute the in situ production of 6 Li in star form- 
ing structures. Indeed, unless CRs are allowed to escape, 
that is, unless the characteristic escape length is significantly 
smaller than 10 g/cm 2 at high redshift the 6 Li abundance 
would greatly exceed the observed value. 

Clearly a definitive result for the production of 6 Li at high 
redshift or at low metallicity will require a more detailed 
model for the ejection and propagation of CRs in the early 
structure forming Universe. The mechanisms described here 
would certainly produce a plateau at low metallicity. The ab- 
solute abundance level is uncertain but can be tied to a set of 
reasonable chosen physical parameters. 
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